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Abstract 

We show that the mean-field time dependent equations in the $ 4 theory can be 
put into a classical non-canonical hamiltonian framework with a Poisson structure 
which is a generalization of the standard Poisson bracket. The Heisenberg invariant 
appears as a structural invariant of the Poisson tensor. 
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1 Introduction 



The purpose of this paper is to give a formulation of the mean-field dynamical equa- 
tions in a bosonic quantum field theory in terms of generalized Poisson brackets. The 
time-dependent Hartree-Bogoliubov (TDHB) equations are dynamical equations for the 
expectations values of one and two field operators. They have been used to study non- 
equilibrium evolution in $ 4 theory |lj . They are quantum dynamical equations. However, 
we will show that they can be put into a classical non-canonical hamiltonian framework 
with a Poisson structure which is a generalization of the standard Poisson bracket. 

Non-canonical Poisson brackets have been found to be usefull in different contexts, 
including plasma and fluid dynamics [0] and condensed matter physics 0. In relativistic 
field theories, it has been shown how to write the field equations in covariant Poisson 
bracket form Q. In non- relativistic quantum theory, a Poisson structure has been recog- 
nized for the time-dependent mean-field dynamics of a system of fermions P or a system 
of bosons P| . A Poisson structure is especially useful to find structural invariants and to 
study stability properties. We will show how this approach can be developped in quantum 
field theory. 

We consider the $ 4 theory. In section 2, we introduce the dynamical variables of 
interest and we recall the TDHB equations which can be written in a compact form 0. 
We define the Poisson tensor in section 3 and we show the Poisson structure of the TDHB 
equations. Another definition of the Poisson tensor from a Lie algebra is considered in 
section 4. In section 5, we recognize the Heisenberg invariant as a structural invariant 
of the Poisson tensor. Finally, in section 6, we write the linearized dynamical equations 
arround the HB static solution. 

2 The mean-field equations in the <P 4 theory : 

We consider a self-interacting scalar field in a Minkowski metric described by the Hamil- 
tonian H = J d d x7i(x) with 

H(x) = l -Ii\x) + \ (V$(f)) 2 + ^U 2 (f) + . (2.1) 

We work in d spatial dimensions. The constants ttlq and b are respectively the bare mass 
and the bare coupling constant. $(x) and IL(x) are the field operators in the Schrodinger 
representation. The initial state of the system is supposed to be given and equal to a 
statistical operator D(to) which may represent a thermal equilibrium or a non-equilibrium 
situation. At zero temperature, it reduces to a projection operator. 
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A time-dependent variational principle [EJ and the choice of gaussian operators for 
the trial space allow to derive time- dependent mean-field equations for the expectation 
values [7[ EJ. A gaussian density-like operator V(t) is completly characterized by a 
vector a(x, t) and a matrix S(x, y, t) defined by : 



ip(x,t) \ 
-iir(x, t) ) 



where 



a(x, t) = 

w(£ y- A - l 2G ^ 1 & ^ ~iT(x, y, t) \ 
\-iT(y,x,t) -2S(x,y,t)J 

<p(x,t) = ^Tr(V(t)<S>(x)) , 
7r(f,t) = -^Tr(P(t)n(x)) , 
G(f,y,t) = J-Tr(V(t)$(xMy)) , 

7i(t) 

s(f,jr,t) = ^Tr(x>(t)n(£)n(y)) , 

T(x,y,t) = -iyTr(P(t)($(x)n(y) + fi($$(f))) , 



2.2 
2.3 

2.4 

2.5 

2.6 

2.7 

2.8 
2.9 



ra(t) = TrV(t) , 

and, for any operator O : = — Tr(T>(t)0). 

The time-dependent Hartree-Bogoliubov (TDHB) equations for the $ 4 theory can 
be written in the following compact form |7| : 

i a = t w , (2.10 

z S = - [(S + r) H (E — t) — (H — t) H (5 + r)] . (2.11 
where r is the 2x2 matrix 

T = (-i I) (2 - 12 

The vector w and the matrix 7i are defined from the variation of 

<H>= —-Tr(V(t)H) = f d d xS(x,t) (2.13 



5 <H>- 



J d d x w(x, t) Sa(x, t)-\j d d xd d y tr (H{x, y, t)6E(y, x, t)) . (2.14 
S(x,t) is the mean-field energy density. For the Hamiltonian fl2.1|) , we have : 

w^x, t) = ^-A + m„ + ^ (/? 2 (£, t) + ~ G{x, x, t)j ip(x, t) (2.15 



w 2 (x, t) = in(x, t) (2-16) 

Unix, V,t) = ~ (-A + m 2 + b - <p 2 (x , t)+ b - G(x, x, t)j 8{x - y) (2.17) 

H 22 {x,y,t) = ^8{x-g) (2.18) 

H 12 (x,y,t) = , H 21 (x,y,t) = (2.19) 

The static HB solution is caracterized by wq = 0, i. e. in the uniform case 

(mo + ~Po + ~Go(z J z)) ¥ >o = (2.20) 
o 2 

and the gap equation 

H = -r coth(/3H T) (2.21) 



which can also be written as 



/JWo = irln(|^) • (2.22) 

More explicitly, in the uniform case, by working in the momentum space, 

coth Pu){p) \ 



with u 2 (p) = p 2 + m 2 + §p§ + §G (f, x). 



In the following we will show that the mean-field equations (|2.10|) and (|2.11|) have 
a Poisson structure. 



3 The Poisson tensor and the Poisson structure of 
the mean-field equations 

We consider first the uniform case. A more general construction of the Poisson tensor 
in the non-uniform case will be given in the next section. We gather in a single vector 
= 1,2,3 the expectation values < 1 >,<x,- and Eij(p), where ^ij(p) is the Fourrier 
transformed of S(x, y, t) : 

y) = J^- d s*(p) • (3.1) 
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We define the Poisson tensor C^ u = {z^, z u } in the following way : 



C 22 - T- • 


/^23 









cflM = 




4 




n(p) 


cgu(p) = 


c 2 3 i 3 h(p) = 




4 


Su(p) 


cll,M = 


8S 12 (p) 








c£u(p) = 


-8H 12 (p) 








^12,22 (p) = 


^21,22 (Pj = 


4 




22 (p) 


^22,12 IP) = 


^22,21 (p) = 




4 


=22 (p) 


CguCp) = 


Cg 2 i(p) = 










(3.2) 



The other matrix elements of this tensor are equal to zero. C 22 is a constant in the 
momentum space, while C 33 depends on the momentum p. This Poisson tensor is a linear 
function of the dynamical variables Eij(p). Let us note that it does not depend on the 
value of the condensat a. 

The Poisson structure associated to the Poisson tensor is defined according to 

for any functions / and g of the dynamical variables z^. 

We check that the operation {/,<?} = & which associates to two functions / and 
g a third function /c is bilinear and antisymmetric. It satisfies also Jacobi's identity and 
Leibnitz's derivation rule ||. (When the Poisson structure is not defined from a Lie 
algebra, we have to check these properties.) 

In particular, we have 

{a „ 9}= i C ^=i n ^. (3.4) 

dOLi 



From the definition Wi = J^, we can write the first TDHB equation ( 2.10Q as : 



Oi = . (3.5) 

We have also : 



= (3.6) 



If we choose g = £, with the definition ( |3.2|) of C 33 and the definition ( 2.14 ) of the matrix 
TC, we obtain : 

i{Su,5} = -4S 12 (p) H 22 (p) 

i{Z 12 ,£} = 2Z 11 (p)H 11 (p)-2Z 22 (p)H 22 (p) • (3.7) 
i{Z 22 ,£} = 4Z 12 (p) H n (p) 
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The second TDHB equation fl2.11| ) can be put in the form : 



(3. 



The TDHB equations in the $ 4 theory can therefore be written as classical dy- 
namical equations, the mean-field energy density £(a, E) playing the role of a classical 
Hamiltonian with the non-symplectic Poisson structure defined by 

It is actually more useful to consider the free-energy density 



S(a,E) - /T 1 ^ 



(3.9) 



as a classical Hamiltonian instead of S. The von-Neuman entropy density S is defined by 

(3.10) 



TrVlogV , _ n 
S = — — h log Tru 



TrV 

For a gaussian density matrix it can be written as a function of the matrix E as : 



5 = / 



(2tt) 3 



!*.{(! (l- T)a (l- T )-l) log(|±I)} 



4 

+- logdet I - ( 
2 5 4 1 -1 



S(1-t) 



1 f 1 

2 ll 



or 



S 



d 3 p 



I V7i OK (^±i) + 1 log(J-i) 



where I is the Heisenberg invariant 



VT-r 



/(ar, y) = J (4 < $(x)$(«) >< II(z)I%) > 

- < $(x)IL(z) + TL(z)3>(x) > < $(z)fi(y) + 5(y)$(z) > 



(3.11) 



(3.12) 



, (3.13) 



y) = J d 3 z (-Su(x, z) E 22 (z, y) + E 12 (x, z) E 12 (z, y)) . (3.14) 
For a pure state, we have : Hn E 22 — E\ 2 + 1 = or I(x, y) = 6(x — y). 



We have {S(E),a} = and {<S(S),S} = 0. Therefore, the TDHB equations 
and (|3.8| ) can be written as : 

di = {ati, F(a,E)} , (3.15) 

Ey = {E ij ,J 7 (a,E)} . (3.16) 

The free-energy density !F(a,E) is a constant of the motion like S(a,E) and S(E). The 
static solution a 0; ( see ec L- (|2.20| ) and ( p.21| )) of the TDHB equations is a minimum of 
T (at zero temperature, it is a minimum of £). In section 6, we will derive the linearized 
mean-field equations arround the static solution. 
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4 Definition of the Poisson tensor from a Lie algebra 



In this section we consider the non uniforme case and we work with the space indices. We 
define the following operators : 

O 1 = 1 (4.1) 

°^ ={-%)) (4 ' 2) 
s r = / 2$(f)$(y) -i($(f)fl(y) + n(y)|.(f))\ 

{x,y> \-z(mfi^) + ^)Hy)) -2n(f)n(y) J l ' j 

where O = O— < O > for any operator O. 

It is easy to show that these operators satisfy a Lie algebra : 

[O", O v \ = f% v O a (4.4) 

More explicitly : 

[0 2 (x),0%(y, z)} = f d d u £(/ 2 23 )™ fe (£, V, z; u) 2 m (u) + (f?) iJk (2, y, z) O 1 (4.6) 

[0%{x, y),OUz, t)] = / <fu <f w E(/! 3 )^(^ £ 3 *1 «, ^ (4-7) 

The structure constants f 22 , /| 3 , /| 3 are delta functions in position space while f 23 depends 
also on the expectation values < 2 (x) >. The other structure constants fg" are equal 
to zero. For instance, we have : 

f 22 (x,y)=r5(x-y) , (yf) 1>12 (f , y, 5) = -2^(f) <5(a? - - 5) . (4.8) 

We then consider the expectation values ^ =< > and we define the Poisson 
bracket C" = {z^ 1 , z v } from the Lie algebra characterized by the structure constant / : 

= {z», z v } = f^z a (4.9) 

Similarly to the structure constant /, the tensor C" has space indices. 

In the uniform case, by working in the momentum space, we can show that the only 
non-vanishing commutators between the operators can be written as : 

Pl0 2 ]=r l3 (4.10) 
PI 0%(p = 0)] = E(/2 3 )™, °l + (fi 3 kjk O 1 (4.11) 

m 

[OUp),OUp)\ = E(/f)^ 3 mn (p) (4.12) 

ran 

where the structure constants / do not depend on the momentum p. In this way we find 
again the Poisson tensor defined in section 3. 
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5 Heisenberg Invariant 



Another presentation of the Poisson tensor can be useful to find the structural invariants 
. We still consider the uniform case and we work in the momentum space. By chosing 
the coordinates < $ >, < II >, < > (p), < till > (p), < $n + fi$ > (p), the Poisson 
tensor can be written as a 5 x 5 matrix, build from two blocks : the symplectic matrix 

and the 3x3 antisymmetric matrix defined by the elements : 



1 
-1 



/?33 



1 .3,22 (P) 



/?33 
u ll,22 



e {< > (p), < $n + n$ > (p)} = 4 < $$ > (p) 

= {< $n + ni> > (p), < nil > (p)} = 4 < nil > (p) 
{< M > (p), < nil > (p)} = 2 < $n + n$ > (p) , 



(5.1) 

(5.2) 
(5.3) 



that is 



M 



( 2 < $5 + ni> > (p) 4 < > (p) \ 

-2 < $5 + ni> > (p) -4 < nn > (p) 

-4 < > (p) 4 < nn > (p) 



(5.4) 



The rank of the Poisson tensor is a function of the dynamical variables. Near a 
regular point, the rank of the Poisson tensor is constant. We can find coordinates (Q,P,C) 
such as the matrix M can put in the form : 



(5.5) 





( 


1 




M' = 


-1 










V 








The Poisson brackets in this coordinates system are given by the symplectic form 

If ] - d/ d 9 df dg 
U,g f dQ dP dP dQ ' 



(5.6) 



We have therefore locally a canonical Hamiltonian system which depends on one param- 
eter C. The function of the dynamical variables which is equal to C is a Casimir function: 
it satisfies {/, C} = for any function / of the dynamical variables 

The coordinate C corresponding to the eigenvalue of the Poisson tensor defined 
by the matrix (|5.4j ) is : 



c(p) = 4 < > (p) < nn > (p)- < $n + n$ > (p) < $n + n$ > (p) . (5.7) 

We recognize the Heisenberg invariant ( 3.13Q and we have : 

C(p) = (5.8) 



9 



6 Linearized equations 



At the HB solution Zq which is a regular point of the Poisson tensor, the matrix defined 
by: 



dz^dz u zo 

is positive definite. The free-energy density T being conserved by the flow defined with 
the Poisson brackets (|3.3|) , this implies that Zq is a stable equilibrium point. 

The TDHB equations write : 

1 8T 

a i = {a i ,T} = -C%{~)-£- (6.2) 



1 „„„ dT 



{E ij ,F} = - i C% kl (E)— (6.3) 



where we have explicitly written the dependence of the Poisson tensor on the matrix E. 
Let us develop these equations arround the static HB solution ao, H at the first order in 
5a = a — ao and 5E = E — E . Because ao, S is a minimum of T , the variation of C 
arround ao, E does not contribute. By introducing 

^\a,E) = ^\a ,~ ) + \ 5z»5z» ^\ Zq , (6.4) 

we obtain : 

^ = {5a t ,^\a,E)}o = ic??(S ) (6.5) 
at i J oybcij) 

where the indice for the Poisson bracket indicates that the Poisson tensor C has to be 
evaluated at the point So- In this way, we obtain the small oscillations equations arround 
the static HB solution also called RPA equations. These equations are of Hamilton type 
with the free energy developped to second order !F^ 2 \a, E) playing the role of the hamil- 
tonian and with C(Eq) as Poisson tensor. !F^ 2 ' is a constant of the linearized equations of 
motion and has a local minimum at ao, Eq. Therefore the linearized dynamical equations 
at Zq have as stable equilibrium. 



The small oscillations equations ( |6.5|) and Q6.6Q can be written as : 
.dSat 



i- 



dt 



^ 2 (Ho) [Bf k 5a k + J Bf kl {q) 8E lk (q)^ (6.7) 



^-^dT 1 = O( s o)(?) ( BfUl) Sam + J ^ B% mn (Q, k) SE nm (k)) (6. 
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The expression of the matrix elements of B^ v in the $ 4 theory is given in appendix. In 
equations (6.7) and Q6.8Q , we have considered small oscillations arround the uniform static 
HB solution : 

a(p, t) = ao5(p) + 8a(p, t) (6.9) 
E(pi,p2,t)=Ea(pi)8(jti+p 3 ) + 8E(pi,p2,t) . (6.10) 

where 

Efafat) = J d 3 xd 3 ye-^-^y~(x,y,t) . (6.11) 

We have then introduced the total momentum P = pi + p 2 and the relative momentum 
Q = — P2) an d we have written the equations for a total momentum P — 0. 

Equations (|6.5|) and ( |6.6|) could also have been obtained by linearizing the original 
expressions ( |2.10|) and (|2.11| ) arround the static solution : 

i8a = t5w (6-12) 

i5E = 2(5E Ht-tH5E) + 2(H 5Ht-t5H S ) . (6.13) 

We define the matrices t,T,r,R (which are the analogs of the RPA matrix TZ of reference 
||) by expanding < H > up to the second order, or the vector w and the matrix TC up 
to the first order : 

8iUi(x, t) = J d d y t i:j (x, y, t) 8atj(y, t) - - J d d yd d z T iJk (x, y, z, t) 8E(z, y, t) (6.14) 

8Hij{x, y,t) = J d d z r ijtk (x, y, z, t)8a k (z, t) - - J d d zd d u R ijtk i{x, y, z, u, t) 8E ik (u, z, t) . 

(6.15) 

We therefore obtain a relation between these matrices and the matrix B of the second 
derivatives of the free energy at the minimum. This relation involves the Poisson tensor 
evaluated at S . 

From the equations fl6.7|) and (|6.8|) , we can study the existence of bound states at 
finite temperature both in the symmetric phase and in the broken phase. This will be a 
generalization to finite temperature of the results obtained by Kerman and Lin at zero 
temperature fTT| . 

To conclude, we have shown that the time- dependent mean-field equations in $ 4 
theory can be written as classical dynamical equations with a non-symplectic Poisson 
structure, the free-energy density playing the role of a classical Hamiltonian. The Heisen- 
berg invariant appears as a structural invariant of the Poisson structure. This Poisson 
structure can be useful to discuss stability properties of trajectories characterized by the 
condensat and the fluctuations of the quantum field at finite temperature. The occurrence 
of a vanishing eigenvalue for the matrix B at a critical temperature for a given coupling 
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constant is the signal of a phase transition from the phase with broken symmetry to the 
symmetric phase. 
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Appendix 

In this appendix we give the expression of the elements of the matrix B in the case 
of the $ 4 theory. 

e;; = m? + ^S + |G„(f,f) , Bg = -1 , Bg = Bg = (A.l) 

^-ott^ (A - 2) 



B 5<='< M > = mJ)L( q ) - \ mJmM s3{p - q) (A ' 3) 

B -* »> = 4^) (sSi + I^bi^O - «> (A ' 6 > 

where x = and 

^ 2 (P) =P 2 + ^o + ^o + ^ (f,f,t) (A.8) 
The other elements of the matrix B are equal to zero. 



References 

[1] O. Eboli, R. Jackiw and S. Y. Pi, Phys. Rev. D37 (1988) 3557 

[2] A. Weinstein, p. 3 in "Fluids and Plasmas : Geometry and dynamics" , J. E. Marsden 
(Ed) Contemporary Mathematics, Vol.28, Amer. Math. Soc, Providence, RI,1984. 

[3] I. E. Dzyaloshinskii and G. E. Volovick, Ann. of Phys. 125 (1980) 67 



12 



[4] J. E. Marsden, R. Montgomery, P. J. Morrison and W. B. Thompson, Ann. of Phys. 
169 (1986) 29 

[5] R. Balian and M. Veneroni, Ann. of Phys. 195 (1989) 324 

[6] M. Benarous, Ann. of Phys. 264 (1998) 1 

[7] C. Martin, Phys. Rev. D52 (1995) 7121 

[8] R. Balian and M. Veneroni, Ann. of Phys. 187 (1988) 29 

[9] P. A. M. Dirac, "Lectures on Quantum Mechanics" , Belfer Graduate School of Science 
Monographs, Yeshiva University, New York, 1964. 

[10] A. K.Kerman and C. Y. Lin, Ann. of Phys. 241 (1995) 185; frep-th/ 9805061 



13 



